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ABSTRACT 

We report the first theoretical estimate of the nuclear-spin dephasing time T<i owing 
to the spin interaction with the two-dimensional electron gas, when the latter is in the 
integer quantum Hall state, in a two-dimensional heterojunction or quantum well at 
low temperature and in large applied magnetic field. We establish that the leading 
mechanism of dephasing is due to the impurity potentials that influence the dynamics 
of the spin via virtual magnetic spin-exciton scattering. Implications of our results for 
implementation of nuclear spins as quantum bits (qubits) for quantum computing are 
discussed. 
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1. Introduction 

Recent ideas [1-3] of utilizing nuclear spins in semiconductor quantum wells and 
heterojunctions as quantum bits (qubits) for quantum computation have generated new 
emphases in the studies of nuclear-spin relaxation and, especially, quantum decoherence, 
in such systems. In this work we consider the case of the integer v = 1 quantum-Hall 
state [4] . The two-dimensional electron gas is then in a nondissipative state. Since the 
electrons mediate the dominant interaction [1,5,6] between nuclear spins, it is reasonable 
to expect that relaxation times of the latter, as well as decoherence/dephasing effects, 
will occur on large time scales. 

Solid-state proposals for quantum computation [1-3] with nuclear spins are all 
presently theoretical. Related proposals to utilize quantum dots [7-15] are also at present 
all in the theory stage. Usually, non-zero-nuclear-spin atoms will be considered placed 
[1] by modern "atomic engineering" techniques in a host material of zero nuclear spin 
isotope. In order to allow positioning with respect to other features of the system, such 
as gate electrodes [2], and making replicas [1], etc., the nuclear-spin separation will be 
larger than the atomic size, typically, of order 20 to lOOA. At these separations, the 
direct magnetic dipole-dipole interaction of the nuclear spins is negligible. 

The dynamics of the nuclear spins is governed by their interactions with each other 
and with their environment. In the regime of interest, these interactions are mediated by 
the two-dimensional electron gas. Various time scales are associated with this dynamics. 
The relaxation time T± is related to energy exchange and thermalization of the spins. 
Quantum mechanical decoherence and dephasing will occur on the time scale T 2 . The 
latter processes correspond to demolition of quantum-mechanical superposition of states 
and erasure of phase information, due to interactions with the environment. Generally, 
there are many dynamical processes in the system, so the times T\ and T 2 may not 
be uniquely, separately defined [16,17]. Theoretically and experimentally, it has been 
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established that processes of energy exchange are slow at low temperatures, so Ti is very 
large, but there still might be some decoherence owing to quantum fluctuations. Gener- 
ally, for various systems, there are extreme examples of theoretical prediction, ranging 
from no decoherence to finite decoherence [18,19,20] at zero temperature, depending on 
the model assumptions. 

In order to consider control ("programming") of a quantum computer, we have to 
identify the time scale T ext of the single-spin rotations owing to the interactions with 
an external NMR magnetic field. We also identify the time scale T int associated with 
evolution owing to the pairwise spin-spin interactions. The preferred relation of the 
time scales is Ti, T 2 T ext , Tj nt , which is obviously required for coherent quantum- 
mechanical dynamics. 

The aim of this work has been to advance theoretical understanding of the time- 
scales of interest for the quantum computer proposal [1] based on nuclear spins in a 
two-dimensional electron gas, with the latter in the integer quantum Hall effect state 
obtained at low temperatures, of order IK, and in high magnetic fields, of several Tesla, 
in two-dimensional semiconductor structures [4]. This system is a promising candidate 
for quantum computing because the nuclear spin relaxation time T\ can be as large as 
10 3 sec. In the summarizing discussion, Section 5, we discuss and compare the values of 
all the relevant time scales. 

Our main result, presented in Sections 2 through 4, is the first theoretical calculation 
of the nuclear-spin dephasing/decoherence time scale T 2 for such systems. We note that 
the recent study [21-23] of the nuclear-spin relaxation time T l5 has relied heavily on 
the accepted theoretical and experimental views of the properties and behavior of the 
electronic state of the two-dimensional electron gas in the quantum Hall regime. These 
electronic properties have been a subject of several studies [4-6,21-29]. We utilize these 
results in our calculation as well. 
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2. The Model 

We consider a single nuclear spin coupled to a two-dimensional electron gas in a 
strong magnetic field, B, along the z axis which is perpendicular to the two-dimensional 
structure. Assuming nuclear spin-|, for simplicity, we write the Hamiltonian as 

H = H n + H e + H ne + H imp (1) 

Here the first term is the nuclear spin interaction with the external magnetic field, 
H n = —^ n Ba z , where 7 n includes h and the nuclear g-factor, and o z is a Pauli 
matrix. 

The second term is the electronic component of the total Hamiltonian (1). Within 
the free-electron nomenclature, the Fermi level lies in between the two Zeeman sub- 
levels of the lowest Landau level. The spin-up sub-level is then completely occupied, so 
the filling factor is v = 1, while the spin-down sub-level is completely empty; note that 
the relevant effective electronic ^-factor in typically negative. In fact, the calculation 
need not be limited to the lowest Landau level. Here, however, to avoid unilluminating 
mathematical complications, we restrict our attention to the lowest level, as has been 
uniformly done in the literature [24-27]. 

The last two terms in (1) correspond to the nuclear-spin electron interactions and 
to the effects of impurities. These will be addressed shortly. The magnetic sub-levels 
are actually broadened by impurities. At low temperatures, the v = 1 system is in 
the quantum Hall state. The interactions of the two-dimensional electron gas with the 
underlying material are not shown in (1). They are accounted for phenomeno logically, 
as described later. 

The electron-electron interactions are treated within an approximate quasiparti- 
cle theory which only retains transition amplitudes between Zeeman sub-levels. The 
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elementary excitations of the electron gas are then well described as magnetic spin ex- 
citons, or spin waves, [24-27]. The spin excitons are quasiparticles arising as a result 
of the interplay between the Coulomb repulsion of the electrons and their exchange 
interaction. A creation operator of a spin exciton with a two dimensional wave vector k 
can be written in terms of the electronic creation operators a* in the spin-down Zeeman 
sub-level and annihilation operators b in the spin-up sub-level as 



L X Ly 



Here £ = y 'ch/eB is the magnetic length, and the p-summation is taken in such a 
way that the wave number subscripts are quantized in multiples of 2n/L y . Note that 
expression (2) assumes the Landau gauge, which is not symmetric under x <-> y. For 
our purposes, the following expression for the dispersion relation of the excitons [24,25] 
provides an adequate approximation, 

. /7T\i/2 fe 2 \£ 2 k 2 
^ = ^(2) U)-2- (3) 

Here A = \g\/j,BB, where /ig is the Bohr magneton, and g is the electronic g-factor, and 
e is the dielectric constant of the material. It has been pointed out [6,23] that the gap A 
in the excitonic spectrum suppresses nuclear spin relaxation at low temperatures. The 
electronic Hamiltonian can be written in terms of the spin exciton operators as 



H e = S + J2 E ^ A l A ^ (4) 
k 

where the c-number So is the spin-independent ground state energy of the electron gas. 
This description of the electronic gas is appropriate only for low density of excitons, 
which is the case in our calculation, as will be seen later. 

- 5 - 



We now turn to the third term in (1), the interaction between the electrons and 
nuclear spins. It can be adequately approximated by the hyperfine Fermi contact term 



H ne = —^ngtlBln • (r e ~ Rn) (5) 

e 

Here hl n and hS e are nuclear and electronic spin operators, respectively, and r e are the 
electron coordinates. The nuclear coordinate R n can be put equal to zero. Such an 
interaction can be split into two parts 



H ne — -ffdiag + -f^offdiag (6) 

where ifdiag corresponds to the coupling of the electrons to the diagonal part of nuclear 
spin operator I n , and -ffoffdiag — to its off-diagonal part. 

The diagonal and off-diagonal contributions can be rewritten in terms of electronic 
creation and annihilation operators as 



(87r/3)7 n ^ B K(0)| 2 x ^ a) / f f \ 

H diag = .1 9 ><r,(a k a q -b k b q ) (7) 



^offdiag = r- L ^ { ' [ a b l a i + ° a L b i ) ( 8 ) 

* y k,q 

Here a ± = \ (a x ±ia y ). The interactions of the electrons of the two-dimensional gas 
with the underlying material are incorporated phenomenologically through the dielectric 
constant and g-factor, see (3), et seq., and via |iuo(0)| 2 and d in (8) above. The latter 
is the transverse dimension of the effectively two-dimensional region (heterojunction, 
quantum well) in which the electrons are confined. The quantity w$(Q) represents 
phenomenologically the enhancement of the amplitude of the electron wave function 
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at the nuclear position owing to the effective potential it experiences as it moves in 
the solid-state material. It is loosely related [5,6] to the zero-momentum lattice Bloch 
wavefunction at the origin. 

For the purposes of the calculations performed here, with the relevant states being 
the ground state and the single-exciton states of the electron gas, one can show that the 
terms in (7) that correspond to different k and q do not contribute, while the remaining 
sum over k becomes a c- number, representing the Knight shift of the polarized electrons. 
Thus ifdiag can be incorporated into the nuclear-spin energy splitting, redefining the 
Hamiltonian of the nuclear spin as H n = \Ta z , where V = 7 n (B + -^Knight)- Note that 
the Knight shift can be used to estimate the value of the phenomenological parameter 
| wo (0)| from experimental data. The off-diagonal coupling (8) can be expressed on terms 
of the excitonic operators (2) as follows [5,6,23], 



The summations over k x cincl ky are taken over all the integer multiples of 2ir/L x and 
27r/L y , respectively. 

The last term in (1) describes the interaction of the electrons with impurities and 
plays a crucial role in nuclear relaxation in the systems of interest. This interaction can 
written in the spin-exciton representation as [23,26] 





where 



C = 



(87r/3) 7n ff/i B | W o(0)| 2 



(10) 



#imp = (2i/L x L y ) U (q) sin V 2 ( k xQy - M*) /2] A^A 



k+q 



(11) 



where U (q) = J U- linp (r) e iq ' r <i 2 r is the Fourier component of the impurity potential 
for electrons in the two-dimensional plane. We will assume [23,26] that the impurity 
potential has zero average and can be modeled by the Gaussian white noise completely 
described by its correlator, (Ui mp (r) U- imp (V)) = Q5^ (r — r'). 

In summary, the relevant terms in the full Hamiltonian (1) can be expressed solely 
in terms of the nuclear-spin operators and spin-excitation operators as 

H = ~T<r x + J2 E k A W + I> k ( A U~ + A * a+ ) + X^^k^k+q (12) 

k k k,q 

where the explicit expressions for £k, (7k and (/>k, q can be read off (3), (9) and (10), 
respectively, and the quantity V was introduced in the text preceding Eq. (9). 

3. Energy Relaxation 

In order to set the stage for the calculation of T 2 , let us first briefly summarize in 
this section aspects of the calculation of the nuclear-spin relaxation time Ti, along the 
lines of [22,23]. The dominant mechanism for both processes at low temperatures is the 
interactions with impurities. Thus, both calculations are effectively zero-temperature, 
single-spin; these assumptions will be further discussed in Section 5. 

We assume that initially, at time t = 0, the nuclear spin is polarized, while the 
excitons are in the ground state, (0)) = |— ) ® |0), where |— ) is the polarized-down 
(excited) state of the nuclear spin and |0) is the ground state of spin-excitons. Since the 
Hamiltonian (12) conserves the ^-component of the total spin in the system, the most 
general wavefunction evolving from |$ (0)) can be written as 

|* (£)> = a (t) |-) <g> |0> + J2 A (t) 1+) ® |lk) (13) 

k 
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with |+) corresponding to the nuclear spin in the ground state and | Ik) describing the 
single-exciton state with the wave vector k. Equations of motion for the coefficients a 
and /3k can be easily derived from the Schrodinger equation: 



iha=^Ta + J2g^ (14) 



2 



q 

In order to solve the system of equations (14)-(15), we introduce Laplace transforms, 
f(S) = J °° f(t)e- st dt, which satisfy 

iSha — ih = -Ta + (7k/3k (16) 

k 

iShpu = ~T0k + #k/?k + ^ </>k, q ^ q + #k« (17) 

q 

Let us first solve (16)-(17) for the case when the interaction of spin-excitons with 
impurities is switched off, i.e., 0k,q = 0. After some algebra we obtain 

1 ( 18 > 



a(s) h ish + r — Ek 

k 

where we have shifted the variable: s = S + iT/(2h), which only introduces an nonin- 
teresting phase factor. 

In the absence of the hyperfine interaction, i.e., for = 0, a(s) in (18) has only 
the pole at s = 0. When the interaction is switched on, the pole shifts from zero. 
This shift can be calculated in a standard way, within the leading order perturbative 
approach, by taking the limit s — > 0, so that ihs ++ r _ Ek — > V r _^ Ek — (r — E^), 
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where V denotes the principal value. This type of approximation is encountered in 
quantum optics [30] . The relaxation rate and the added phase shift of the nuclear-spin 
excited-state probability amplitude a(t) are given by the real and imaginary parts of 

2 

the pole, TjT- = ^ Xlkfi'k^ (r — -E'k) and uj = r-g k res P ec tively, so that oc 

e -t/(2Ti)+iwt j t j g Q^viQug that d ue to the large gap in the spin-exciton spectrum 
(3), r A, the energy conservation required by the delta function above can never be 
satisfied, and so in the absence of interaction with impurities, Ti = oo. It also transpires 
that Ti is infinite [30], as will become apparent later. 

Interactions with impurities, described by the last term in (12), will modify the 
solution of (16)-(17), and, as a consequence, the energy conservation condition. In par- 
ticular, if the impurity potential is strong enough, it can provide additional energy to 
spin-excitons, so that their energy can fluctuate on the scale of order T thus making 
nuclear-spin relaxation possible. This mechanism [22,23] corresponds to large fluctua- 
tions of the impurity potential U (r), which usually occur with a rather small probability, 
so Ti is very large for such systems. 

In order to carry out the above program quantitatively, one has to solve the system 
of equations (16)-(17) with nonzero 0k,q- Such a solution is only possible within an 
approximation. One can introduce the effective spin-exciton self-energy Ek in (18), so 
that ihs+ Y_ E ifts+r-g k +s k • integral equation for Ek can then be derived, taking 
the continuum limit in (16)-(17). Solving this equation would allow one to calculate 
the relaxation rate from (18). However, in order to satisfy the energy conservation, 
we require V — E^ + Ek = 0, so the self-energy should be rather large, of order E^. 
Therefore, as a result of the spectral gap of the excitons, the perturbative approach is 
inadequate as it automatically assumes that |Ek| <C l-^kl- Instead, a certain variational 
approach [23,31] has been adapted to evaluate T l5 consistent with the experimental 
values [33,34] of order 10 3 sec; for further discussed see Section 5. 
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4. Dephasing Mechanism 

We argue that in order to calculate the phase shift due to the impurity potential, 
one can indeed use the perturbative solution of (16)-(17). Indeed, phase shifts result in 
virtual processes that do not require energy conservation and therefore are dominated by 
relatively small fluctuations of the impurity potential simply because large fluctuations 
are very rare. Moreover, the terms in the sum in (18) that contribute to the relaxation 
rate, do not contribute to the phase shift, see the discussion above. This consideration 
also applies when the self-energy is introduced. 

One can show that the contribution to dephasing linear in </>k,p vanishes due to 
symmetry. Thus, let us solve (16)-(17) perturbatively up to the second order in (/>k p 
and perform the inverse Laplace transform of a(s). Within this approximation, the 
pole of a(s) in the complex- s plane is imaginary, so that \ce(t)\ = 1. We conclude that 
a(t) oc e luJut and f3k(t) = 0, where the part of the phase-shift responsible for dephasing 
is 



The zeroth-order term in (19) was dropped as irrelevant for our calculation of the 
dephasing time. Since V is much smaller than i?k, it was also omitted. 

As expected, the perturbative solution does not describe the energy relaxation 
(7\), but it does yield the additional phase shift due to the impurity potential. We 
will see shortly that this phase shift, when averaged over configurations of the impurity 
potential, produces a finite dephasing time, T2. 

Let us consider the reduced density matrix of the nuclear spin, given by 




(19) 



p n (t)=[Tt e \9(t))(9(t)\] u (20) 

recall (13). Here the trace is partial, taken over the states of the spin-excitons, while 
the outer brackets denote averaging over the impurity potential. The trace over the 
spin-excitons can be carried out straightforwardly because within the leading-order per- 
turbative approximation used here they remain in the ground state; all excitations are 
virtual and contribute only to the phase shift. The diagonal elements of p n (t) are not 
influenced by virtual excitations and remain constant. 

The off-diagonal elements of p n (t) contain the factors e ±tuJut . It is the averaging of 
these quantities over the white-noise impurity potential U (r) that yields dephasing of 
the nuclear spin. In order to proceed, let us rewrite (19) more explicitly. From (9)-(ll), 
after changing summation index (k — > k — a ^ £ ) in the first sum in (19) we obtain 



(21) 

Here we use the following shorthand notation for the z component of a vector product 
Q]^ = k x q y — k y q x . 

It is appropriate to assume that impurity potentials are short-range, i.e., a <C £, 
where a is the scale of variation of £/i mp (r). This assumption and the white-noise 
property of the impurity potentials, are required to make the problem amenable to 
analytical calculation. Thus, the main contribution to the Fourier transform U(p), 
dominating the summation in (21), comes from large wavevectors p (and q), of order 
a -1 ^> t~ x . Therefore one can replace the exponent e~ V(p +c i) 2 by the Kronecker symbol 
5 q _ p , to obtain a simplified expression for uu 



Now the sum over k can be carried out because for p ^> k, we can assume that -Ek+p — 
E p ~ E c (£ 2 /2)p 2 , where E c = (tt/2) 1/2 [e 2 /(e£)}. 

Moreover for large p, the factor sin 2 {(£ 2 /2) [k, p] z } can be replaced by its average, 
1/2. Finally, we get 



1 _e-^ k2 sin 2 f [k,p] 2 1 1 r 



k ^k^k+p ~~£ c £ 2 p 2 2W (A + £ c ^ 2 ) 

The integral can be evaluated explicitly; specifically, for <C 1 we get 



-lit 2 
e s K 



fceta ^ ~ — (24) 

o (A + E c ^k 2 ) 2 PE C A 



so that 



2C 2 



E 



[/(p)[/(-p) 



M C 2 £ 4 A (L x L y ) 2 p- 



(25) 



Recall that we have assumed the white-noise distribution for the impurity poten- 
tial, (Ui mp (r) [/i m p (r')) = Q5( 2 ) (r — r'). This corresponds to the following probability 
distribution functional for the Fourier-transformed potential, 



P[U(p)] = iVexp 



^QL X Ly 



^[/(p)[/(-p) 



(26) 



The latter expression, and other approximations assumed earlier, allow to reduce the 
averaging of e lU}ut to a product of Gaussian integrations. The off-diagonal elements of 
the nuclear-pin density matrix are, thus, 



Poi 



IT 



L X LyP 2 



— exp 



■EM 1 



IT 



L X LyP 2 



(27) 



where r = 4QC 2 t/(7i7r£ c 2 rA) 



We are interested in real part of the sum in (27), which represents decoherence/de- 
phasing of the nuclear spin. The off-diagonal elements decay exponentially as 



p 01 ~ exp 



_-ii>( i+ cdkO 

The summation over p in (32) can be converted into integration 



(28) 



5> ( 1 + 7f-r 



L ' L " ' d 2 p\n(\ ■ 



(L x L y ) 2 p*J (2n) 2 J " V (L x L y ) 2 p± 
Explicit calculation then yields the result that p i ~ e _T / 16 or p 01 ~ e~*/ T2 , where 



(29) 



2M 2 A . , 

T - = w (30 » 

with E/a = Q/(2n£ 2 E 2 ). 

5. Results and Discussion 

The quantity C/2 characterizes the strength of the impurity potential with respect 
to the Coulomb interactions [26] . Let us summarize typical parameter values [26] for a 
GaAs heterojuction, which is the system best studied in the literature. For magnetic 
field value B = 10 T, we have the following values of parameters, I = 0.8 x 10 _8 m, 
E c = 3x 10 -21 J, C = 2.5 x 10 -36 Jm, A = 4.6 x 10 -23 J. From experimental data 
for electronic mobility, one then estimates [26] U 2 ^ 0.0025, yielding T 2 ~ 40 sec. We 
emphasize that this is an order of magnitude estimate only, because of the uncertainty 
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in various parameter values assumed and the fact that the parameters, especially the 
strength of the disorder, may vary significantly from sample to sample. For instance, 
there is another estimate of the disorder strength Q available in the literature [23], 
obtained by fitting the value of T\ to the experimentally measured 10 3 sec [33,34], as 
cited earlier. This yields an estimate for T 2 that is smaller, T 2 ~ 0.5 sec. Generally, we 
expect that with typical-quality samples T 2 may be a fraciton of a second or somewhat 
larger. 

Let us point out that to date there are no direct experimental probes of dephasing 
by the disorder-dominated mechanism identified here for dilute nuclear-spin positioning 
appropriate for quantum computing. Such systems where never engineered. For those 
materials whose atoms have nonzero nuclear-spin-isotope nuclei, specifically, for GaAs 
(spins 3/2), we are aware only of one experiment [35] where indirect information on 
dephasing can be obtained from the linewidth. However, in that case the dipolar in- 
teractions cannot be neglected and likely provide the dominant dephasing mechanism. 
For quantum computing, the host material will have to be isotope-engineered with zero 
nuclear spins, e.g., Si [2]. 

Let us now compare various time scales relevant for quantum computing applica- 
tions. The relaxation time T\ is of order 10 3 sec [23,33,34]. For the spin-spin interaction 
time scale T; nt , values as short as 10 -11 sec have been proposed [1,5]. These estimates 
are definitely overly optimistic and require further work. Since such calculations re- 
quire considerations beyond the single-spin interactions, they are outside the scope of 
the present work. For T ext , modern experiments have used NMR field intensities corre- 
sponding to the spin-flip times of 10~ 5 sec. This can be reduced to 10~ 7 sec, and with 
substantial experimental effort, perhaps even shorter times, the main limitation being 
heating up of the sample by the radiation. 

Thus, the present information on the relevant time scales does not show violation 
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of the condition Ti,T 2 ^> T ex t,T in t, stated in the introduction, required for quantum 
computing. To firmly establish the feasibility of quantum computing, reliable theoretical 
evaluation of Ti nt is needed, as well as experimental realizations of few-qubit systems 
engineered with nuclear spins positioned as separations of order 30 to 100 A. 

We also note that typical lab samples, for which the parameter values used were 
estimated, have been prepared to observe the quantum-Hall-effect plateaus in the re- 
sistance. The latter requires a finite density of impurities. However, for the quantum- 
computer applications, a much cleaner sample would suffice. Indeed, as suggested by 
our calculations, T 2 is mostly due to dephasing owing to virtual spin-exciton scattering 
from impurities. Therefore, the value of T 2 can be increased by using cleaner samples. 

We acknowledge helpful discussions with Drs. S. E. Barrett, M. L. Glasser and 
R. Mani. This research was supported by the National Security Agency (NSA) and 
Advanced Research and Development Activity (ARDA) under Army Research Office 
(ARO) contract number DAAD 19-99-1-0342. 
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